Spin-polarized Shapiro steps and spin-precession-assisted multiple Andreev reflection 
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We investigate the charge and spin transport of a vohage-biased superconducting point contact 
coupled to a nanomagnet. The magnetization of the nanomagnet is assumed to precess with the 
Larmor frequency, ljl, due to ferromagnetic resonance. The interplay between the ac Josephson 
current and the magnetization dynamics leads to spin-polarized Shapiro steps at voltages \V\ — 
hjJL/'ien for n = 1,2, ... and the subharmonic steps with n > 1 are a consequence of multiple Andreev 
reflection (MAR). Moreover, the spin-precession-assisted MAR generates quasiparticle scattering 
amplitudes that, due to interference, lead to current-voltage characteristics of the dc charge and 
spin currents with subharmonic gap structures displaying an even-odd effect. 

PACS numbers: 75.76.-|-j, 74.50.-|-r, 75.50.Xx, 75.78.-n 
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Introduction. Control of the electron's spin degree of 
freedom has led to a number of spintronics applications. 
[H These applications often rely on using the sponta- 
neously broken spin-rotation symmetry of a ferromagnet 
to act as a reference for the spin of itinerant electrons. 
In contrast to ferromagnets, the electron-electron inter- 
action in a BCS superconductor leads to a spinless order 
parameter. While these two states of matter usually are 
incompatible in bulk materials, they can be combined 
in nanoscale junctions to create new phenomena. 0^ 
The interplay between superconductivity and ferromag- 
netism in such hybrid superconducting junctions allows 
for probing the spin degree of freedom as well as for co- 
herently controlling spin-polarized currents. In ad- 
dition, superconducting junctions coupled with magneti- 
zation dynamics have been studied [l6l - [l8j . The lowest- 



magnetic impurity 23[ or a single-molecule magnet [3C 
34| . We show that the interplay between the Josephson 



energy magnetization excitation can be accessed by the 
application of an external magnetic field, which starts a 
precession of the magnetization around the direction of 
the field This ferromagnetic resonance (FMR) 

mode can also be achieved by a coupling between the 
magnetization and the Josephson effect: the ac Joseph- 
son current generates an oscillating magnetic field that 
resonantly excites the magnetization precession 23l- 25|. 



In Ref . [24] , the coupling between a spatially dependent 
order parameter and the magnetization dynamics was 
studied in a tunnel junction whose width, W , was larger 
than the superconducting coherence length, ^o, and this 
coupling was detected as a rectification of the ac Joseph- 
son charge current. However, in the regime W < S^q, 
which is the case for a superconducting quantum point 
contact (SQPC) [H, 27[, the spin- rotation symmetry of 
the superconducting order parameter may be broken lo- 
cally and one may anticipate that spin-dependent modifi- 
cations of the superconducting correlations will influence 
the transport properties t2S,i29j. In this Letter, we study 
the coupling of the ac Josephson effect and magnetization 
dynamics of an SQPC containing a nanomagnet, e.g. a 



and magnetization oscillations create Shapiro-like reso- 
nances that lead to a rectified spin current. In addition, 
we find a rich subgap structure of the current-voltage 
characteristics of the dc charge and spin currents dis- 
playing features related to the emission or absorption of 
energy quanta corresponding to the precession frequency. 
However, the features associated with an odd number of 
Andreev reflections are suppressed due to interference. 

The transport properties of nonmagnetic SQPCs can 
be understood in terms of Andreev reflection. An ap- 
plied bias voltage leads to the occurrence of multiple 
MAR [35|] and the ac Josephson effect, which is char- 
acterized by the Josephson frequency, wj — 2eV/h. The 
ac current of SQPCs also includes higher harmonics of 
the Josephson frequency, i.e. j{t) = X^raJn^™"'^* 
ting h= I) [36l-38|. In addition, MAR processes produce 
a subharmonic gap structure of the dc current that can 
be seen in the current- voltage characteristics as steps at 
voltages V = 2A/en, n = 1,2,... The ac Josephson ef- 
fect can be probed by the application of microwave ra- 
diation as their interplay leads to Shapiro steps and a 
rich subgap structure due to photon- assisted MAR [sil - 
4l| . In the case of SQPCs with magnetic interfaces, the 
transport properties are modified due to spin-dependent 
effects 0, E.g., a magnetic interface may produce 

different transmission probabilities for the spin-up and 
spin-down bands, and Pi , which may lead to the oc- 



currence of spin currents [ill IT^ 



Here, we consider a point contact consisting of two su- 
perconducting leads coupled over a nanomagnet. The left 
(L) and right (R) leads are assumed to consist of s-wave 
superconductors in the clean limit (slj. We further as- 
sume that the magnetization of the nanomagnet can be 
treated classically within the macrospin model and hence 
be described by single spin vector, S Under FMR 
conditions, external dc and rf fields are applied to create 
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an efTective magnetic field, H, in which the spin precesses 
with the Larmor frequency = 7|-ff|, where 7 is the 
gyromagnetic ratio. This effective field also includes any 
effects of the resonant coupling between the ac Joseph- 
son current and the magnetization dynamics. The nano- 
magnet's spin, whose dynamics can be described by the 
Landau-Lifshitz-Gilbert equation [i^], generates a time- 
dependent exchange field that is felt by the quasiparticles 
tunneling across the junction and can be incorporated 
into a phenomenological tunnel Hamiltonian that allows 
for spin fiip scattering J,— i>t (t~^i) which is accompanied 
by absorption (emission) of an energy quantum a;^ by 
the quasiparticle (see Ref. [l^ for details). 

The transport characteristics depend on the time- 
dependent superconducting phase difference across the 
contact, ip{t) = ipQ+LUjt, where ipo is the initial phase dif- 
ference, in addition to the Larmor frequency. The spin S 
traces out an in-plane angle described by x(0 = Xo +i^Lt 
in analogy with ip(t). The scattering resulting from 
the MARs combined with the precession-induced spin 
fiips transfers quasiparticles into sidebands with ener- 
gies = e + nuij + muiL, where n = 0, ±1, ±2, . . . and 
m = 0, ±1. The charge (/z = 0) and spin (^ = x,y,z- 
component) currents in lead a = L,R can then be writ- 
ten in terms of all harmonics as 
ijp 



(a) D^=0.8 



(b) D^=0.2 
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Model. The currents are obtained using nonequilib- 
rium Green's function techniques and the calculation is 
based on the Hamiltonian approach in which the 
leads are described by Keldysh Green's functions in the 
quasiclassical approximation [4^ and the interface con- 
taining the nanomagnet is treated as a strong impurity 
44| . The effects of (p{t) are accounted for in the stan- 
dard way by a transformation of the unperturbed lead 
Green's functions, g^, where """ denotes a matrix in 
Keldysh space, that is expressed as Ua = eKjp[iipa(t)T3/2] 
in Nambu space (" ""). The voltage-biased lead Green's 
functions, ga{t,t'), are then expressed in terms of the 
equilibrium ones, ga{t — t'), as g^{t,t') = h{l{t)g^ {t — 
t')Ua{t') where X refers to one of the retarded (i?), ad- 
vanced {A) or Keldysh [K) components. The equilib- 
rium Green's functions have a simple Fourier transfor- 
mation given by ga{t — t') = J{de/2TT) e''''^^^^*"> gais). 
Since the leads consist of s-wave superconductors, the 
retarded component of gds) is given by g^ = T3g^{e) + 
Tj'^is) iay, where g'^je ) = -ne/n, /^(e) = -nA/n, 
H. — -^/l Ap — (e + ?0)^. The advanced component is 
9a{^) — ^i[9a{^)V ^'i while the Keldysh Green's function 
is given by g^(e) ^ [g^{e) - g^{e)] tanh(e/2T). 

The magnetic interface is characterized by a tunnel 
Hamiltonian, Ht = J^ka-k'a' c\ f,„VL,ka\R,k' a' CR^k' a' + 
/i.e., in which the hopping amplitudes, ?)i,i?(i) = VRLit) — 
v{t) = v''{t), depend on the instantaneous direction of the 
spin, S(t). Within the quasiclassical approximation, the 
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FIG. 1. (Color online) Charge background current and 
conductance, Gb = djg/dV, for a vr junction with (a) 
T>i; = 0.8 and (b) T)^ = 0.2. The current is normalized 
by the normal conductance given by the two spin channels, 
Gn ~ [e^/h][V^ + ©J.]. The conductance curves in panel 
(a) have been offset with 3 {Vi = 0.8), 2 {T>i = 0.5) and 
1 (X>^ — 0.2). In panel (b), the conductance is plotted for 
Vf = 0.2, = 0.1 and ul/A = 0.2. In all plots, i? = 7r/8 
and the temperature is T = 0. 





FIG. 2. (a) First-order tunneling process in which a spin-up 
electron-like quasiparticle, e^, absorbs a precession quantum 
cjL and undergoes a spin flip, (b) Second-order process which 
leads to destructive interference of the combined Andreev re- 
flection amplitudes of the hole-like quasiparticles ht and h4,. 



Fermi surface limits of the hopping amplitudes can be 
expressed as Vk;k' = {vo +vs{S{t) ■ (T))5{k — k'). We can 
define the spin-up (spin-down) transmission probability 

as 2?t(i) = 4i^t(i)/[l + ^ta)]^' "^^^^^ ^t(W = VQ±vscosd 
and d is the angle between S and H . 

The hopping elements' influence on the lead Green's 
functions is captured by the t-matrix equation ia{t, t') = 
f„(i,i') + JdtiJdt2tait,h)g^{ti,t2)L{t2,t'), where 

fL/_R(t,i') = VLR/RL{t)9R/L{t,t')vRL/LR{t')- The 

time dependence of the superconducting phase of the 
Green's functions can be transferred to the hopping 
elements by the transformation Ua{t) according to 



UL/R{t)vLR/RL{t)Ul,j^{t) = exp[±i(^(i)f3/2]{)(t). The t- 
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matrix equation then becomes an algebraic equation in 
energy space and 

ia(t t') — g-i(rtyo+mxo)+i(rmj+inaJi,)f' _ ^'^J™ 

{Ut)]n= J ^e--*[i„«r)]. (2) 

By only considering currents on the left side of the junc- 
tion (thereby dropping the index a = L), one can write 
the Fourier components as = [i(£™)] and the current 



where we have defined i^" ^ vq + vsctz cosd, — 
^[cr^c ± icTj,] sin?? and = if > 2. The advanced 
and Keldysh matrices are similarly defined. 

Results. - The dc charge and spin currents of Eq. ([1]) 
consist primarily of contributions from the (j^)o com- 
ponents and can be described as background currents 
in the spirit of Ref. [s^l- There is also the possibility 
of Shapiro-like resonances between the Josephson and 
the Larmor frequencies at voltages = —{m/n)uJL/'2e, 
m ^ 0. The Larmor precession produces a spin structure 
of the t matrices in which the index m = 1(— 1) cor- 
responds to a spin-raising (spin-lowering) matrix. This 
spin structure causes the Shapiro resonances to give con- 
tributions to the spin-polarized current only and thus 
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FIG. 3. (Color online) The z component of the dc background 
spin current. In all plots, i? = 7r/8 and T = 0. 



components of Eq. ([T]) as 

= I |^Tr{«:^[Cg° - 5™C]<}, (3) 

where = 5(0' ^° = efa, = diag(crj, cryCTi(T.y)/2 
for i — x,y, z. The i-matrix components can be found 
from the equation C = f ;7 + + K'n+A+i + 

B^'^_ii\^_i} where the retarded components are given 

by' 



(4) 
(5) 
(6) 
(7) 

I 

result in a Shapiro spin current given by Jghapiro = 
En,m=±i e-'("^"+'"'^''Hj')n '5(V"- V„"). This bchavior is 
in contrast to that of a weak link coupled to a nanomag- 
net whose magnetization performs a precessional motion. 
In that case, the resulting time-dependent magnetic field 
produces Shapiro- type resonances that give contributions 
to the charge current [isj . 

Figure [1] shows the current-voltage characteristics as 
well as the conductance for the charge background cur- 
rent, jg. As can be seen, the most prominent feature 
is the subharmonic gap structure appearing at voltages 
eV — (2A ± WL)/n, where n = 1,3,... This structure 
can clearly be seen in the conductance where it is dis- 
played as distinct satellite peaks for junctions in the tun- 
nel limit. The subharmonic gap structure is similar to 
that of the current- voltage characteristics generated by 
photon-assisted MAR which shows features due to the 
absorption or emission of photons [39| . However, in the 
case of MAR processes influenced by Larmor precession, 
a quasiparticle can only absorb or emit one quantum of 
ljl since these processes are accompanied by a spin flip 
(see e.g. Fig. [2] for the absorption processes with n — 1 
and 2). As a result, there is only one satellite peak on 
each side of the feature corresponding to a tunnel pro- 
cess of order n. In addition, only the processes for which 
n is an odd number display side peaks. One can show 
that this suppression is due to destructive interference 
of the Andreev reflection amplitudes corresponding to 



i=-2 

2 



and [BX,;1, = E e~ 

j=-2 



[g Jn-l 
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the side peaks of processes with n = 2,4,... Defining 
SV^ = [el^— 2A±wl]/2A, the height of the current steps 
at 6V^ ^ 1 can be approximated in the tunnel hmit by 
j| « e ((5F±) eAf sin^ -dl {[eV ± ujl]/2A), where 9 is 
the Heaviside step function, I (a) = [2aE{^yl — 1/a?) — 
K{^1 — l/a2)/a] and K and E are the complete elliptic 
integrals of the first and second kinds. In the limit <C 
2A, one obtains / ([eF ± wl]/2A) « f [l + ^SV^] and 
the heights of the current steps •d,tV= {2A ± ujl)/^ are 
thus oc v1 sin^ §. 

The z component of the background spin current, j'^, 
is shown in Fig. [3] (the in-plane components are zero). 
For transmission probabilities T>\ ^ T>,i and = 0, our 



results reproduce those of Ref. jl2| for zero spin mix- 
ing. In the case of nonzero Larmor frequency, the spin 
current is finite even for = V^, as can be seen in 
panel (b). In the tunnel limit, the spin current can be di- 
vided into a spin-filter current and a spin-pump current, 
Jb = j'mtcr + Jpump- The z component of the spin-filter 
current is given^ = (A/27r)(Pt - V^)I{eV/2A) 

for eV > 2 A [12|. The spin-pump current can corre- 
spondingly be approximated by jp^jjp^^i = ±j±/2e. For 
eV > 2A + LUL and ujl <C A, the total spin-pump current 
is given by jpu,„p^^ = {3/4:)v^ sin^ cjl if wl ^ A and is 
hence oc and does not depend on the bias voltage. 

Now, we turn to the Shapiro-like resonances in the spin 
current. Since the Larmor precession requires m to take 
the values m — ±1 for the Shapiro resonance condition, 
Jshapiro ^P'^^ polarized in the xy plane (see Ref. [2^ 
for a detailed analysis of the t matrices' spin structure). 
Defining {r)nSiV -V„) ^ E™=±i[(j"):r + (r)-„]^(^- 
T^'") where Vn — ±Wi/2en, one can write the nth com- 
ponent as (jshapiro)n = R^j'')nR S{V -Vn) , whcrc n > 1 
and the transformation R = exp[^(— tt-^sq i Xo)H ■ cr] ro- 
tates the vector (j'')„ around the direction of the external 
magnetic field, H, through an angle —nipQ -I- sign(V^i)xo- 
The existence of a dc current that is spin polarized in the 
xy plane implies that the rotational symmetry around 
the z axis is broken and that the magnitude of the 
Shapiro currents depend on the initial angle of nanomag- 
net's magnetization direction. This behavior is analo- 
gous to the dependence on the initial superconducting 
phase difference which also is present in the behavior 
of microwave-irradiated Josephson junctions [39]. The 
Shapiro contributions strongly depend on the transmis- 
sion probability as can be seen in Fig. |4l 

The iJL values in Fig. [T]are chosen for clarity but may 
be smaller in practice. Consider an SQCP consisting of 
Nb, whose superconducting gap is A ^ 1 meV, contain- 
ing a nanomagnet that under FMR conditions reaches 
precession frequencies oi uj^ ~ 10 GHz which is well be- 
low the critical magnetic field. He, of Nb ^ISj- Then, 
typically wl/A = 0.01 and the temperature is restricted 
to T < 100 mK. Increasing H closer to the critical field 
value decreases A and hence increases the ratio ujl/A 
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FIG. 4. (Color online) The Shapiro-resonance contributions 
(a) (j»)i, (b) {r)2, and (c) {fh for = = V. Panels 

(d) -(f) show the Shapiro-resonance contributions (d) {j")!, 

(e) (j'')2, and (f) {j'')s for = 0.5. In all plots, i9 = tt/8 
and r = 0. 



which allows for a better resolution of the subgap struc- 
tures oi jg and jg. Alternatively, using the ac Joseph- 
son current to resonantly excite S{t), lol values corre- 
sponding to an effective magnetic field \H\ > He can be 
achieved f23l-[25j . Detection of the spin-polarized Shapiro 
currents would then be a measurement of the FMR fre- 
quency of the nanomagnet. 

Conclusions. - In conclusion, we have calculated the 
dc charge and spin currents through a voltage-biased 
superconducting point contact coupled to the spin of a 
nanomagnet that under ferromagnetic resonance condi- 
tions precesses with the Larmor frequency uj^. We have 
shown that coherent multiple Andreev reflection in the 
presence of the Larmor dynamics leads to a subharmonic 
gap structure of the dc charge and spin currents with 
features at eV = (2A ± uji^)/n originating from the ab- 
sorption or emission of a precession quantum. However, 
destructive interference suppresses these features for even 
values of n. In addition, the Larmor dynamics combined 
with broken spin-rotation symmetry generates Shapiro- 
like resonances resulting in an additional spin current 
polarized in the xy plane. 
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